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ABSTRACT 

Galaxies over four decades in stellar mass are shown to obey a strong correlation between u-K colors and 
atomic-gas-to-stellar mass ratios (G/S), using stellar mass-to-light ratios derived from optical colors. The correla- 
tion holds for G/S ranging from nearly 10:1 to 1:100 for a sample obtained by merging the SDSS DR2, 2MASS, 
and HYPERLEDA Hi catalogs. This result implies that u-K colors can be calibrated to provide "photometric 
gas fractions" for statistical applications. Here this technique is applied to a sample of ^35,000 SDSS-2MASS 
galaxies to examine the relationship of gas fractions to observed bimodalities in galaxy properties as a function of 
color and stellar mass. The recently identified transition in galaxy properties at stellar masses ~ 2-3 x 10 ll) M 
corresponds to a shift in gas richness, dividing low-mass late-type galaxies with G/S ~ 1 : 1 from high-mass galax- 
ies with intermediate-to-low G/S. Early-type galaxies below the transition mass also show elevated G/S, consistent 
with formation scenarios involving mergers of low-mass gas-rich systems and/or cold-mode gas accretion. 

Subject headings: galaxies: evolution 



1. INTRODUCTION 

Analyses of galaxies in the Sloan Digital Sky Survey (SDSS) 
have demonstrated two distinct bimodalities in galaxy proper- 
ties: a bimodality between recent-burst dominated and more 
continuous star formation histories (SFHs) as a function of stel- 
lar mass M», divided at ~ 3 x 10 10 M Q (Kauffmann et al. 
2003b), and a bimodality between blue late-type and red early- 
type galaxy sequences as a function of optical color, divided at 
u-r ~ 2.2 (Strateva et al. 2001; Hogg et al. 2002; Blanton et al. 
2003b). Recently, Baldry et al. (2004) have partially unified 
these observations, demonstrating a color transition within each 
of the two galaxy sequences atM* ~2x 10 10 M , as well as an 
increase in the relative number density of red sequence galax- 
ies above ^2-5 xlO 10 M Q . They also argue that the number 
density of the red sequence is consistent with a major-merger 
origin. However, the cause of the color and SFH transitions at 
~2-3 x 10 10 M Q remains to be explained. 

Several physical processes that influence SFHs may imprint a 
transition mass on the galaxy population. Supernova-driven gas 
blow-away will preferentially affect halos with small escape ve- 
locities (Dekel & Silk 1986), although simulations suggest that 
the baryonic mass threshold for blow-away may be closer to 
10 7 M Q than to 10 ll) M (Mac Low & Ferrara 1999). Cold- 
mode gas accretion may dominate in low-mass halos whose gas 
fails to shock to the virial temperature (Birnboim & Dekel 2003; 
Katz et al. 2003); here analytic estimates give a threshold mass 
of a few times 10 11 M including dark matter, so a link to the 
observed transition at M* ~ 2-3 x 10 10 M is plausible. Finally, 
observations suggest that inefficient star formation may be typ- 
ical of disk-dominated galaxies with V c < 120km s" 1 , possibly 
reflecting the relative importance of supernova feedback as op- 
posed to other turbulence drivers in supporting the interstellar 
medium against gravitational instability (Dalcanton et al. 2004). 

All of these processes involve gas - its expulsion, accretion, 
or rate of consumption. Thus examining how the gas properties 
of galaxies vary with color and stellar mass may offer vital clues 
to the origin of the transition mass and the color shifts within the 
red and blue sequences. Unfortunately, tracing the dominant 
neutral phase of the interstellar medium requires H I 2 1 -cm line 
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observations, which are challenging even at the modest redshifts 
probed by the SDSS. To make full use of the statistical power 
of the SDSS, an alternate strategy is required. 

Building on earlier optical work (e.g., Roberts 1969), Bothun 
(1984) has shown a remarkably tight correlation between Hi 
mass-to-H-band luminosity ratios and B-H colors. Going one 
step further, the present work describes a method for estimat- 
ing atomic-gas-to-stellar mass ratios using u-K colors from the 
SDSS and Two Micron All Sky Survey (2MASS) databases. 
This "photometric gas fraction" technique is calibrated using 
H I data from the recently expanded HYPERLEDA H I catalog. 
When the technique is applied to a sample of ^35,000 SDSS- 
2MASS galaxies at z < 0.1, the transition mass of 2-3 xlO 10 
M is observed to correspond to a shift in gas richness found 
separately in both galaxy color sequences. This result implies 
that any explanation of the transition mass via gas physics must 
directly or indirectly affect both early- and late-type galaxies. 

2. DATA 

Optical, near-IR, and H I data were obtained from the SDSS 
second data release (DR2, Abazajian et al. 2004), the 2MASS 
all-sky extended source catalog (XSC, Jarrett et al. 2000), 
and the HYPERLEDA homogenized H I catalog (Paturel et al. 
2003). Merged catalogs were constructed containing all z < 0.1, 
r < 17 .77, K < 15 galaxies with positions matched to within 6" 
and with reliable redshifts and magnitudes based on data flags 
and cataloged errors (magnitude errors < 0.3 in K, < 0.4 in H I, 
and < 0.15 in ugr). The 2MASS magnitude limit was set fainter 
than the completeness limit to improve statistics on dwarf and 
low surface brightness galaxies. As the 2MASS XSC has un- 
even depth, it probes significantly fainter than the complete- 
ness limit in some areas of the sky. Because of their marginal 
detectability, galaxies with H 1-derived gas-to-stellar mass ra- 
tios greater than two were targeted for individual inspection, 
and eight were rejected as having unreliable 2MASS or SDSS 
pipeline reductions. These rejections exacerbate the shortage of 
IR-faint galaxies. The final samples are: SDSS-HYPERLEDA 
(575 galaxies), SDSS-2MASS-HYPERLEDA (346 galaxies), 
and SDSS-2MASS (35,166 galaxies). An additional require- 
ment for the SDSS-2MASS sample was that the Local Group 
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motion-corrected redshift be greater than 1000 km s . 

All optical and IR magnitudes used here are fitted magni- 
tudes, i.e. SDSS model magnitudes and 2MASS extrapolated to- 
tal magnitudes. The SDSS magnitudes are corrected for Galac- 
tic extinction using the DR2 tabulated values and k-corrected to 
redshift zero using kcorrect v3.2 (Blanton et al. 2003a), while 
the 2MASS ZT-band magnitudes are k-corrected using k(z) = 
-2.1z (Bell et al. 2003). Distances are computed in the concor- 
dance cosmology f2 m = 0.3, Qa = 0.7, Hq = 70 km s" 1 Mpc" 1 . 

3. RESULTS 

Fig. la shows the basic correlation between M-band and 
21-cm apparent magnitudes m u and mm for the SDSS- 
HYPERLEDA sample. Its existence is not surprising: M-band 
light is a tracer of young massive stars, and the birth rate of 
young stars is known to depend on the available gas reservoir 
(as in the global correlation between disk-averaged star forma- 
tion rate and gas surface density, Kennicutt 1989). The presence 
of young massive stars may also enhance Hi detection (e.g., 
Shaya & Federman 1987). The absolute magnitude correla- 
tion obtained by distance-correcting m u and m H i is of course far 
stronger than the correlation in Fig. la, but at the cost of non- 
independent axes. In any case, what is relevant for predicting 
m H i from m u is not correlation strength but scatter. Most of the 
0.92 mag scatter in the m„-m H i relation is not explained by the 
errors. This scatter likely represents variations in M-band extinc- 
tion, molecular-to-atomic gas ratios, and the physical conditions 
required to convert a gas reservoir into young stars. Even with- 
out calibrating these factors, the m„-m H i relation is sufficiently 
tight for the present application. 

Figs, lb and lc plot atomic-gas-to-stellar mass ratios 
(G/S) against u-r and u-K colors for the SDSS-2MASS- 
HYPERLEDA sample. Gas masses are derived from H I fluxes 
with a helium correction factor of 1.4, and stellar masses are 
derived from Zf-band fluxes using stellar mass-to-light (M/L) 
ratios estimated from g — r colors as in Bell et al. (2003). The 
resulting correlations are distance-independent and extremely 
strong, with Spearman rank correlation coefficients of 0.75 and 
0.69 for u — K and u-r respectively. Note that the calibration 
sample spans the color-M* relation (Fig. 2a), and the color-G/S 
relation is tighter than the color-M* relation for these galaxies 
by ^25% in both u-K and u-r. The strength of the u-K 
color-G/S relation derives both from the underlying m„-m H i re- 
lation and from the close correspondence between Zf-band light 
and stellar mass. The latter correspondence is assumed within 
this work and may not apply to all starbursting systems (Perez- 
Gonzalez et al. 2003); however, Kauffmann et al. (2003a) find 
that spectroscopically determined M/L ratios generally agree 
well with color-based M/L ratios, even in the low-mass regime 
where starbursts are common. 

The large dynamic range of the u-K color-G/S relation 
makes the relation forgiving of errors and thus well-suited to 
low-precision estimation of photometric gas fractions. The 0.37 
dex scatter in the relation provides a basis for error estimation. 
Furthermore, galaxies of low and high mass define broadly sim- 
ilar u-K color-G/S relations in their regime of overlap (tri- 
angles and dots in Fig. lc). It should be borne in mind that 
the generality of the photometric gas fraction technique as cur- 
rently formulated relies on the fact that heavily dust-enshrouded 
star formation, as in luminous infrared galaxies, is rare in the 
low-z universe (Sanders & Mirabel 1996). In dusty systems 
one might find high G/S linked to red u-K colors. Also, the 



calibration given here could significantly underestimate actual 
gas-to-stellar mass ratios if stellar M/L ratios are much lower 
than assumed and/or molecular gas corrections are large (a con- 
troversial topic; see Casoli et al. 1998; Boselli et al. 2002). The 
M/L ratios used here are roughly consistent with maximum disk 
assumptions for spiral galaxies (Bell et al. 2003). Within the 
ugriz magnitude set, the best alternative to the u-K color-G/S 
relation is the u-r color-G/S relation. Its larger scatter may 
in part reflect the fact that TT-band magnitude errors will move 
points along the u-K color-G/S relation, but away from the 
u-r color-G/S relation. However, the effect from cataloged er- 
rors is quite small (arrows in Figs, lb and lc), so the greater 
u-r scatter seems to be mostly physical. 

Fig. 2a plots u-r color vs. for the — 35,000-galaxy SDSS- 
2MASS sample, with points color-coded to indicate photomet- 
ric gas fractions G/S pri0 t (computed from u-K colors and the 
fit in Fig. lc). The well-known red and blue sequences (refer- 
ring to galaxy color, not color-coding) are roughly separated by 
the dashed line (an approximation to the separator of Baldry 
et al. 2004). Within each sequence, the color-coding reveals 
a shift in gas fractions near a threshold mass of Ml <~ 1-3 x 
10 10 M Q . Massive red-sequence galaxies are extremely gas- 
poor (G/S p hot as low as 1:100, red/yellow color-coded points), 
whereas for red-sequence galaxies below M[, intermediate gas 
fractions (G/S pri0 t ~ 1:10, green points) are the norm. Likewise, 
massive blue-sequence galaxies have intermediate gas fractions 
(G/S p hot ~ 1:10, green points), but blue-sequence galaxies be- 
low M\ are typically gas rich (G/S p /, 0( ~ 1:1, blue points). 

These results are shown in binned form in Fig. 2b, where the 
contours show the conditional probability distribution of u-r 
on M*. In this plot, a vertical slice through the contours at a 
given mass gives the one-dimensional probability distribution 
for u-r colors at that mass. The conditional probability distribu- 
tion is formed by weighting each galaxy by 1/V mat , where V max 
is the maximum volume within which it could have been de- 
tected, then normalizing the counts in each color x mass bin by 
the total counts in each mass column. Bins with fewer than four 
galaxies are not considered. This algorithm is most robustly ap- 
plied to truly magnitude-limited samples, but the results shown 
here look qualitatively similar to those of Baldry et al. (2004) 
for a magnitude-limited SDSS sample. Each galaxy is assigned 
the smallest of three V max estimates, based on (i) the magnitude 
limit in r, (ii) the magnitude limit in K, and (iii) the distance 
limit imposed by the z < 0. 1 selection requirement. 

The color-coded symbols in Fig. 2b illustrate gas fraction and 
morphology trends using l/V mflA -weighted bin averages. Ovals 
and S-shapes identify early and late type morphologies, respec- 
tively, based on the concentration index C r (defined as rgo/r^o, 
where r^o and rsa are the 90% and 50% Petrosian radii). A 
value of C r = 2.6 is commonly used to divide early and late types 
when true morphological type information is lacking (Strateva 
et al. 2001; Kauffmann et al. 2003a; Bell et al. 2003). Fig. 2b 
indicates transitional values (2.55 <C r < 2.65) with both type 
symbols. Some of the bluest low-mass galaxy bins show tran- 
sitional types, which will not be examined here. Otherwise, the 
blue sequence forms a broad, low C r strip, with a shift in gas 
richness at Ml ~4-3 x 10 10 M Q . At a similar mass, the red 
sequence shows a coordinated shift in gas richness and C r . 

Intriguingly, the increased gas content of the red sequence 
below Ml causes the red and blue sequences to fuse (nearly) in 
a plot of GISphot vs. M* (Fig. 2c). The double-peaked, slightly 
S-shaped contours of the G/S p hot vs. M* distribution thus look 
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FIG. 1 . — (a) Apparent n-band magnitude vs. apparent H I magnitude for the SDSS-HYPERLEDA sample. A bisector fit yields m u = 0.33 + 0. 98mm (thick line) 
with cr = 0.92 mag (thin lines), (b and c) Atomic-gas-to-stellar mass ratio G/S vs. u-r and u-K color for the SDSS-2MASS-HYPERLEDA sample. Bisector fits 
yield log(G/S) = 1.46— 1.06(n — r) and log(G/5) = 1.87 — 0.56(w — K) (thick lines) with cr = 0.42 dex and 0.37 dex, respectively (thin lines). Arrows indicate the 
effect of a 0.3 mag error in K (the maximum allowed by selection). The color bar in panel c provides a key to the photometric gas fractions in Fig. 2. Horizontal 
lines demarcate gas-rich, intermediate, and gas-poor regimes. Dots and triangles mark galaxies with M* above and below 2 X 10 Mq, respectively. 
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FIG. 2. — Galaxy color vs. stellar mass for the SDSS-2MASS sample. Vertical lines show the transition mass interval at 1-3 X 10 10 Mq. (a) Individual u — r 
data points color-coded by photometric gas fraction G/S p ; lo , (computed using u-K). The dashed line divides the red and blue sequences (referring to galaxy color, 
not color-coding), approximating the optimal divider of Baldry et al. (2004). The natural separation of the sequences is clearest in the high-resolution online figure. 
Black dots show the Hi calibration sample of Figs, lb and c. (b) Gas fraction and morphology trends within the red and blue sequences. Contours show the 
conditional probability distribution of u—r on M„, in bins of 0.25 X 0.5 in color xlog(M„ ). Each galaxy is inverse- volume-weighted (see text) and all bins of a given 
mass are normalized by the total for that mass column. Each bin has three slightly offset symbols whose colors show the weighted mean and ±lcr values of G/S ; ,/, 0( , 
with errors propagated from the scatter in the u — K color-G/S relation. Ovals and S-shapes indicate early and late type bins based on weighted mean concentration 
index (see text), (c) G/S p /, , (or equivalently u—K) vs. stellar mass. Small points show individual values, while contours give the conditional probability distribution 
of G/S pkot aIM1 u-K on M„, in bins of 0.5 X 0.5 in color xlog(M„). Large points show weighted mean G/S p ; lo , trends for the red and blue sequences separately, 
using the u-r division in panel a. All contours start at 0.04 and increase by 10 X at each step. 



qualitatively similar to the bimodal distribution of SFH vs. M* 
reported by Kauffmann et al. (2003b). 2 This similarity sug- 
gests that the bimodality in SFHs may be intimately related to 
changes in G/S. The contours in Fig. 2c show a broad shift near 
~10 10 Mq, while the individual trends in the red and blue se- 
quences change slope over the range ~l-3 x 10 10 Mq (though 
the latter changes are not independent of the method of separat- 
ing the two sequences). Kauffmann et al. find a single-sequence 
bimodal structure in C r vs. M* as well, perhaps because of the 
decreasing C r of the red sequence below M\. 

Despite low C r values and high gas contents, the low-mass 
red sequence does fit into the general merger scenario for the 
origin of early-type galaxies, in a way that may illuminate the 



Kauffmann et al. results. The low-mass red sequence maps 
closely to an abundant population of faint, moderately gas-rich 
SO and SO/a galaxies observed in the Nearby Field Galaxy Sur- 
vey, a survey designed to represent the natural distribution of 
galaxy types over a wide range of luminosities (Jansen et al. 
2000; Kannappan et al. 2002, see Fig. 3 of the latter). Based 
on the frequency of gas-stellar counterrotation in this popula- 
tion and the scarcity of dwarf intermediate-type spiral galax- 
ies, Kannappan & Fabricant (2001) argue that >50% of low- 
luminosity SO's may form via late-type dwarf mergers. Such 
gas-rich mergers would naturally produce remnants with mod- 
est C r , explaining why low-luminosity early types are predom- 
inantly SO rather than E galaxies. Moreover, low-luminosity 
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SO's quite often have blue, starbursting centers, despite red outer 
disks (Tully et al. 1996; Jansen et al. 2000), and such color gra- 
dients correlate with evidence of recent interactions (in all mor- 
phological types, Kannappan et al. 2004). The SFH measures 
adopted by Kauffmann et al. (2003b), which are based on 3"- 
aperture spectroscopy, may emphasize the starbursting centers 
of this class of galaxies (and conversely, the quiescent bulges 
of many high-mass late-type systems), reinforcing the single- 
sequence structure of the resulting SFH vs. plots. 

4. DISCUSSION 

A connection between star formation histories and gas frac- 
tions is in some sense obvious: gas must be consumed to form 
stars, so old red stellar populations will tend to be associated 
with diminished gas supplies. However, galaxies also accrete 
and expel gas, so this simple view misses much of the story. 
A complete picture must explain why gas fractions depend on 
mass, and in particular why there is such a close coincidence be- 
tween the transition to gas richness (atomic-gas-to-stellar mass 
ratio <~ 1:1 in the blue sequence) and the shift to recent-burst 
dominated SFHs below ~ 2-3 x 10 10 M Q . Critical transition 
masses are predicted by scenarios involving starburst-driven gas 
blow-away, inefficient star formation below a gravitational in- 
stability threshold, and/or cold-mode gas accretion (Dekel & 
Silk 1986; Verde et al. 2002; Birnboim & Dekel 2003; Katz 
et al. 2003; Dalcanton et al. 2004). The abundance of gas in 
low-mass (10 s — 10 10 M Q ) galaxies seen in this paper is hard to 
explain if global gas blow-away is a dominant process in this 
mass regime. However, localized gas blowout or strong feed- 
back could inhibit efficient widespread star formation in low- 
mass disk galaxies. Such scenarios would not explicitly ac- 
count for the gas in low-mass red-sequence galaxies, but for- 
mation via mergers of low-mass late-type systems, as discussed 
above, could allow low-mass red-sequence galaxies to acquire 
their modest gas excesses from gas-rich progenitors and thereby 



inherit the progenitors' threshold mass for increased gas frac- 
tions. Alternatively, it is possible that the transition in star for- 
mation modes at ^2-3 x 10 10 M© is not a cause, but an effect of 
changing gas fractions, as in cold-mode accretion scenarios. If 
so, low-mass galaxies may form stars reasonably efficiently and 
still appear gas rich. The excess gas in low-mass red-sequence 
galaxies could in this case represent post-merger cold-mode ac- 
cretion, possibly as part of a process of disk regrowth. 

In conclusion, this paper has demonstrated a shift in galaxy 
gas mass fractions at M\ ~ 1-3 x 10 10 M©, likely related to the 
shift in SFHs observed near the same stellar mass by Kauffmann 
et al. (2003b). The link may be causal in either direction, de- 
pending on the relative importance of supernova blowout, feed- 
back, and cold-mode accretion processes in determining Ml . To 
establish this result, a technique has been introduced to estimate 
photometric gas fractions based on the correlation between u-K 
colors and atomic-gas-to-stellar mass ratios. This correlation is 
interesting in its own right (see also Bothun 1984) and will be 
further examined and applied in future work. 
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